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Abstract: Radio Cherenkov emission underlies detection of high energy 
particles via a signal growing like the particle energy-squared. Cosmic ray- 
induced electromagnetic showers are a primary application. While many 
studies have treated the phenomenon approximately, none have attempted 
to incorporate all the physical scales involved in problems with time- or 
spatially-evolving charges. We find it is possible to decompose the calcu- 
lated fields into the product of a form factor, characterizing a moving charge 
distribution, multiplying a general integral which depends on the charge evo- 
lution. In circumstances of interest for cosmic ray physics, the resulting 
expressions can be evaluated explicitly in terms of a few parameters obtain- 
able from shower codes. The classic issues of Fraunhofer and Fresnel zones 
play a crucial role in the coherence. 

1 Introduction 

Coherent radio Cherenkov emission is a remarkably effective method for de- 
tecting high energy particles. The history of the effect goes back to Jelly, 
who first asked whether cosmic ray air showers might produce a radio signal 
|I[ . Askaryan subsequently predicted a net charge imbalance in air show- 
ers, and coherent radio power scaling like the energy of the shower squared. 
Substantial radio emission from atmospheric electromagnetic cascades was 
observed more than 30 years ago. 0, [IJ. Progress in ultra-high energy air 
showers has sparked renewed interest, and new observations of radio pulses 
have been reported recently ||. The current pilot project RICE [|, || uses 
radio Cherenkov emission to detect lOOTeV and higher energy neutrinos in 
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Antarctic ice. The radio Cherenkov signal is the most efficient known mech- 
anism H [17|, H 0, H for detecting neutrinos of lOOTeV and above in solid 
media, yielding detection volumes of order lkm 3 per radio detector for PeV 
neutrinos on ice targets. At PeV energies and above, the neutrino interaction 
cross sections offer fascinating new tests of Standard Model physics || and 
new physics [10]. Tomography of the Earth is also possible with Pe\^-scale 
neutrinos [|ll|]. Radio Cherenkov signals have also been used to search recently 
for neutrinos and cosmic rays with energies upwards of 10 20 eV impinging on 
the MoonQ. 

Cherenkov radiation is also an intrinsically interesting and beautiful phys- 
ical phenomenon. Coherence is a basic feature of electrodynamics, and the 
coherent enhancement of Cherenkov radiation in the microwave region has 
been observed in the laboratory [14|]. The Argonne wake- field accelera- 



tion project |TJ] has successfully generated extremely large microwave field- 
strengths by manipulating coherent radiation from an intense electron beam. 

Despite a long history, the previous literature apparently does not con- 
tain a careful treatment of evolving charge distributions, such as those of 
electromagnetic showers in air or ice, which incorporates all important fea- 
tures of the problem. The problem is intricate because of a multitude of 
scales. When an electromagnetic shower evolves, it produces a pancake of 
charge with a finite thickness, a finite width, probed at a finite wavelength 
of radiation, and for a finite distance over which the shower is big. All this 
occurs at a finite distance from the detector. Results on evolving and finite- 



sized charge distributions are few. Tamm |15| grappled with the problem 
of a charged particle on a track of limited length in the early days of the 
theory. Askaryan anticipated a coherence cut-off in air showers at high 
frequencies of order the inverse pancake size, imposed somewhat by hand. 
Allan gave physical arguments and order of magnitude estimates based on 
one of Feynman's electrodynamic formulas. Kahn and Lerche |16|] attempted 
to resolve the coherence issue using superpositions of infinite tracks. McKay 



and Ralston |30| and Alvarez-Muniz and Zas {AZ) [18j considered the in- 



fluence of the LPM effect at ultra-high energies. Zas, Halsen, and Stanev 
(ZHS) ||19|| , and AZ reported results from summing asymptotic far-fields 
track-by-track in Monte Carlo calculations of great complexity. 

We present an approach which incorporates all the scales and allows a 
general analysis. Main results include an expression for the electric field in 
a factorized form. The "factorization" occurs when distance scales can be 
separated: the characteristic size of the moving charge distribution must be 
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substantially smaller than the scale over which the charge develops. This 
condition is well satisfied for all cosmic ray applications we have examined. 
A form factor characterizes the moving charge distribution, which multiplies 
a charge-evolution integral. Not all the scales decouple: subtleties coming 
under the classic description of Fraunhofer and Fresnel zones need careful 
treatment. Finally, the generic situation can be summarized by analytic 
formulae. This is indispensable given the large parameter space. For exam- 
ple, the numerous and varied numerical plots obtained from immense Monte 
Carlos can be summarized by a few parameters. With the parameters fixed, 
predictions can be made for any number of circumstances. 



2 Consequences of Coherence 

Before beginning analysis we review a few basics. The well-known Frank- 
Tamm (1937) formula uses an exact solution to Maxwell's equations for a 
uniformly moving charge on an infinitely long track. The solution is kine- 
matic and can be obtained by the trick of boosting the charge to a speed 
faster than light in the medium. Extension to a track of finite length has 



pitfalls. Tamm's 1939 finite-track formula fig] assumes a uniform charge e 
traveling at a uniform velocity v along the z-axis for —L/2 < z < L/2. 
Tamm gives the energy loss dP per angular frequency duj per solid angle dfl: 

d 2 P ne 2 2 . 2n sin 2 X 

(ujL) 2 sm 2 9——, (1) 



dudn 4tt 2 c 3V ' X 2 

where X = nuL/(2c)(cos6 c — cos 9) and cos# c = c/(nv); n is the index of 
refraction.^] This formula has been cited in the high energy physics litera- 
ture, and used to interpret experiments observing Cherenkov radiation in the 



millimeter wavelength range [13 



Tamm's finite-track formula includes two competing and distinct phys- 
ical processes: the Cherenkov radiation of a uniformly moving charge, and 
bremmstrahlung or acceleration radiation from charges modeled as starting 
and stopping instantly at the track's endpoints. The interference of the sud- 
den start and stop contributions with the straight line contribution leads 
to strong oscillations in the angular distribution. Compared to a typical 
high energy process the acceleration at the endpoints is fake, that is, the 



1 We diverge from current practice and use the special symbol c to represent the speed 
of light, otherwise known as "1". 
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Tamm model is unreliable. This is because a charge created by pair pro- 
duction is accompanied by an opposite charge, which coherently shields the 
pair from radiating until the oppositely charged partners gradually separate. 
This "Perkins effect" has been observed |2(J and is closely related to the 
coherence phenomena in QCD of color transparency ||21| . At the end of the 
Cherenkov processes, charges also do not stop instantly, but instead slow 
gradually to subluminal speeds. While the slowing has stochastic elements, 
it is better approximated by a uniform deceleration than by a catastrophic 
disappearance of charge. Of course, the evolution of a cosmic ray shower over 
many radiation lengths and involving billions of particles is an even smoother 
macroscopic process. Thus the Tamm formula and related approximations 
may misrepresent the physics if the artificial treatment of the endpoints play 
a major role. Conversely, experimental situations with conditions close to 
those assumed by the Tamm formula can be constructed: Takahashi et al. 
Ul3|| report on the sudden appearance of a charge in a cavity with metallic 
boundary conditions, leading to a strong mixing of Cherenkov and boundary- 
condition effects. 

2.1 Fraunhofer versus Fresnel 

Our study uncovers another, deeper problem with certain asymptotic as- 
sumptions of the Tamm-type approach. In a typical application of Cherenkov 
radiation in high energy physics, we might have a track of length L ~ lm, 
observed at a distance R ~ lm, and in the optical regime with R/ A ~ L/X > 
10 6 . The application to radio detection in ice might have R ~ 100 — 1000m, 
A ~ 0.1 — lm, L ~ 10m, with L/X ~ 10 — 100 with R/X even greater. Cosmic 
ray air showers develop and are observed over tens of kilometers. In all cases, 
all lengths are large in units of the wavelength. Given L large enough for 
the acceleration contributions to be small, the Tamm formula, Eq. (0), might 
appear ideal at first sight. Indeed for ujL — > oo, the sin 2 X/X 2 distribution 
approaches 27rc/(nuL)5(cos6 — cos6 c ). Integrating over angles we recover 
the well-tested Frank- Tamm result for an infinitely long track, 

d 2 P _^\ in 2 Qc 
dudL c 2 

Yet the Tamm formula is quite inapplicable to such problems. This is 
evident from the formula's prediction that the radiated energy will be con- 
centrated in coordinate space at 9 = 9 C , up to a small width due to diffraction. 
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For many of the physical situations cited, the energy is actually spread rather 
uniformly over the length of a cylinder surrounding the charge's trajectory. 
This is a broad angular distribution extending over angles A6 ~ L/R, where 
R is the distance to the receiver. The Tamm formula, or any asymptotic far 
field approximation, does not depend on the distance R, and cannot describe 
this simple truth. True enough, the momenta (wave numbers) of photons 
have directions that may be peaked at 9^ ~ C , but this is not the same thing 
as the power density dP/dfl seen on a sphere surrounding the system going 
like S(cos8 — cos# c ). There is no paradox: if one fixes L arbitrarily large, 
and then moves to an asymptotically distant location R — > oo, the photons 
traveling at the Cherenkov angle will appear to come from a point source 
whose angular size is diffraction-limited. The Tamm formula is derived by 
taking the limit R — > oo; once taken, the case of finite L/R is unavailable. 

To see this breakdown from a different perspective, one can use simple 
dimensional analysis and geometrical reasoning. The Fourier transform of 
the electric field E u has dimensions of mass. Cherenkov radiation for the 
"long uniformly moving track" has cylindrical symmetry. At cylindrical ra- 
dius p from the track, the energy per length 2irpE 2 (p) remains constant. 
Dimensional analysis plus the cylindrical symmetry forces the electric field 
to go like I p. This is nothing like the usual radiation from accelerated 
charges which has fields falling like 1/R in three-dimensional space. And this 
peculiarity applies out to arbitrarily large distance p, provided the track is 
long enough. But if p is taken so large that the radiation appears to emerge 
from a point source, the E field must fall like Lu/p ~ Luj/R. (The factor of 
L comes from the linear power per unit length dependence. Momentarily we 
will examine this in more detail.) 

The breakdown of Tamm's formula is thus due to an interchange of limits. 
Tamm's formula is obtained by making the Fraunhofer approximation, which 
fails under a broad range of finite track-lengths. The Fresnel zone describes 
a complementary far-field region where the Fraunhofer approximation must 
be modified. The basic physics of the Fresnel zone for Cherenkov radiation 
is elementary but requires some care. 

2.2 The Coherence Zone 

Consider (Figure 1) a charge moving on a straight line. Let R(t) be the 
instantaneous distance from the charge to the observation point. Information 
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propagating in the medium at speed c m will arrive simultaneously from the 
track if dR/dt = c m . This is the Cherenkov condition: dR/dt = vcosO = c m 
for velocity v oriented at angle 9 relative to the direction R. Note that R(t) 
is the radius from the charge to the point, not the vector position. 

Due to the geometry of the track and observation point, uniform motion 
produces acceleration of R(t). If dR/dt were constant, the fields arriving 
would all be in phase for the whole track length. However, the acceleration 
relative to the observation point produces an extra radial change of order 
AR = l/2(d 2 R/dt 2 )(At) 2 . Coherence of modes of wavelength A is then 
maintained only over a finite region of AR < A. Since d 2 R/dt 2 = v 2 sin 2 9/R, 
we solve to find the condition At coh < \J R\j [v sin 6). Equivalently, there is 
a finite spatial coherence region for given R, given wavenumber k = 2tt/\, 
namely 

Az coh < ^R/(ksm 2 6) 

over which the "sonic boom" of radiation is built coherently. Since Az co h ~ 
v^R, the coherence zone grows to infinite size as R — > oo: but this limit 
cannot be taken carelessly. 



2.3 Coherence of Evolving Charge Distributions 

We now return to the emission from an electromagnetic shower or other time- 
evolving charge distribution. To a reasonable approximation, the number of 
particles in a highly relativistic shower scales like the primary energy divided 
by a suitable low energy threshold. The charge imbalance near the shower 
maximum is of order 20% e of the total number of particles. These numbers 
have been confirmed over and over, with each generation of numerical simula- 
tion contributing further detail. The origin of the emitted Cherenkov power 
going like the shower-energy squared is basic electrodynamics: the electric 
field will scale like the charge, and the radiated power scales like the electric 
field- squared. 

The evolving shower has a finite length scale a over which it is near its 
maximum, and radiating copiously. This length scale, known as the "longi- 
tudinal spread" in cosmic ray physics f22 |, is akin to the length scale L of 



Tamm's approach but represents a smooth onset and decline of maximum 
power. The shower maximum-length scale a is determined by the material, 
and is conceptually distinct from the shower's total depth to reach the max- 
imum (which goes like the logarithm of the energy) or the charged pancake 
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size (which is fairly constant once the shower is developed.) A cartoon of 
these ideas is given in Figure 1. 

There are then two characteristic limits. Suppose the longitudinal spread 
of the shower is "short" compared to the coherence length, a <C Az co h- Then 
coherence is maintained over the whole range where the current is appre- 
ciable. The amplitude is proportional to the total length a over which the 
current was "on", times 1/R. This is then the R — » oo limit, or Fraunhofer 
approximation, and looks like normal radiation. From dimensional analysis 
(and here we recall the discussion earlier), E w ~ auo/R. 

However, in the limit a > Az co h, the coherence length is not as big as the 
longitudinal spread, and coherence only exists over the smaller of the two. 
Adding amplitudes only over the region Az coh and weighted by 1/R, the E u 



field goes like uAzch/R = Ju/R, This behavior is rather different from the 
previous case: indeed Cherenkov radiation is fundamentally a Fresnel-zone 
effect, as seen by the l/\fR dependence of the fields. 

Both the Fresnel and Fraunhofer limits are far-field approximations, in 
the sense that kR ^> 1 is assumed. The subtlety lies in the dimensionless 
ratio 

kn 2 

T] = (a/Azcoh) = — sin 9 
ri 

which controls how the limit R —>■ oo is taken. Confusion on this point is 
easy; one has R/X ~ 10 4 in the same regime, and yet R is not large enough 
for a "large R" Fraunhofer approximation to apply, exactly because the term 
"large R" is undefined until the limit parameter rj is specified. In the RICE 
experiment one typically has a ~ 1 — 2m, u ~ 100 — 1000MHz, and R ~ 
10 3 m, so 7] < 1 holds. Extension to closer observation points, or to energies 
where the LPM effect can give a much larger a, makes 77 ^> 1 possible. Partly 
due to the obscurity of the coherence criteria, the Fraunhofer approximation 
has received much attention in the previous literature [|1^, except for 



those estimates using fields with "cylindrical" symmetry [23, 17, 24, 25 



2.4 General Set-Up 

Let E^x) be the time- Fourier transform of the E field, with similar notation 
for other fields. The Maxwell equations for a dielectric medium are 
V • I) w — 47rp, cV x B w = AnJu — iuD^, V • B w = 0, cV x E w = iuB w , where 
D^x) = e^E^ix). There is a wave equation for the vector potential A%(x), 
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given by c(V 2 + k 2 )A^(x) = —4ttJ^(x), with k = ujyfe/c. Then we have 

cA w (x) = fd 3 x' 6XP f dt ' exp {iut')J{t\ x'). (2) 

*J 00 00 *J 

The 4-potential = (A , A) has been defined in a generalized Lorentz gauge 
cV-A+edA° jdt = appropriate to the medium. The 4-current J M = (p/e, J), 
p. Since the components of J are related by J = vp, we have A = A°ev/c. 
We calculate A w and then use A^ = v ■ Ac/(ev 2 ). For radiation problems 
the denominator factor l/\x — x'\ is replaced by 1/R. This is standard, with 
corrections of order a 2 /R 2 or similar effects in the "near field" regime, which 
is not our subject. 

2.5 Evading The Fraunhofer Approximation 

The Fraunhofer approximation is the textbook expansion for the phase 
exp (ik\x — x ; |) w exp (ik\x\ — ikx-x'), dropping terms of order k\x'\ 2 /R. 
All existing simulations make this approximation for the phase, and for good 
reasons: the subsequent integrations become much simpler. The integrand in 
Eq. (Q) oscillates wildly. A Monte Carlo simulation has to find the surviving 
phases from myriad cancellations, due to the phases generated over the length 
of each track, and then summed over thousands to millions of tracks moving 
in three dimensions. For a ITeV shower the code of ZHS runs in about 
20 minutes on a workstation. Increasing the energy by a factor of 100, the 
calculational time scales up faster than linear, and computer time becomes 
prohibitive. For this reason various strategies to rescale the output have 
been used in arriving at the published values of electric fields. Even for TeV 
energies, standard Monte Carlo routines such as GEANT challenge a work- 
station's capacity. For cosmic rays of the highest energies the entire approach 
of direct numerical evaluation is unfeasible. 

Unfortunately the Fraunhofer approximation also neglects terms in the 
phases, namely k\x'\ 2 /R, that may be of order unity given our previous 
discussion of length and frequency scales. We must avoid this step. Progress 
is possible due to the translational features of the macroscopic current 
J M (t',x'). A rather general model is 

J(t', x') = vn{z')f{z' - vt', p'). (3) 

An even more general situation will be discussed shortly. The charge packet 
travels with the speed v , chosen here to be along the £-axis of the coordinate 
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system. The function f(z' — vt f , p') represents a normalized charge density of 
the traveling packet, with p' the transverse cylindrical coordinate relative to 
the velocity axis. We normalize / by Jdz'd 2 p' f(z' , p') = 1. In ice the packet 
is about Az' = 10cm thick in the longitudinal direction, and Ap' = 10cm in 
radius (the Moliere radius) in the vicinity of the shower maximum. These 
size scales are limited because of relativistic propagation. The time evolution 
of these scales is negligible near the shower maximum, and indeed the Moliere 
radius is usually approximated by a material constant for the whole shower. 
Similarly, in air showers the scale of charge separation is small compared to 
the scale of shower longitudinal spread. 

The shower's net charge evolution appears in the factor n(z'). With our 
normalization, n(z') represents the total charge crossing a plane at z' . The 
symbol n max will denote the maxiumum value of n(z'); later we will see that 
the electric field scales linearly with n max . The longitudinal spread a is a 
property of n(z') near the shower maximum. The model neglects charge 
(current) left behind, and moving at less than light-speed in the medium, 
which does not emit Cherenkov radiation. We do not have a sharp cut-off at 
the beginning or end of tracks, and the function n(z') will vary smoothly. 

2.6 Factorization for the Fresnel Zone 

Now while we cannot expand around x' — (the Fraunhofer approximation), 
the conditions of the problem do permit an expansion around p ' — (the 
shower axis), namely for R(z') = [(z — z') 2 + p 2 ] 1 ^ 2 , that 



For typical values in this problem, the second term is ~ 10 times smaller 
than the first, and the third is ~ 10 3 smaller than the second. For the 
exponent in Eq. (fj), the third term does not contribute if kAp' 2 /R <C 1, 
that is uo < 250GHz. 

Collecting terms, we have 



x — x 



[(,- 2 ') 2 + (p-pOT /2 



X 



// 




(4) 
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We have shifted the t' integral which produces the translational phase in the 
z' integral. This gives the factorization: 



~ F(q)A FF (r]), 



(5) 



where 



F(q) = jc 



iq-x' 



(6) 



vcRA; 



TFF 



vl 



FF 



(7) 



and 



I FF ( V ,9)= Jdz>eM<t>(z% 



<fi(z') = ik(z'cos9 c + R{z',p)) + \ogn(z'). 



(8) 



Here q = (ou/v,q±), q± = kp/R, and x' = (z',p'). Provided F(uS) <C 1 in 
either frequency region kAp' 2 /R > 1 or kAz' 2 /R ^> 1, the decoupling of 
the integrals is excellent. Here Ap' and Az' refer to the regions over which 
the charge exists near the maximum. F(q) is the form factor of the charge 
distribution, which happens to be defined, just as in the rest of physics, in 
terms of the Fourier transform of the snapshot of the distribution. From our 
definitions F(0) = 1. 

It is worth noting that the dependence on orientation of q is observable. 
For example, in a Giant Air Shower, where the mechanism of charge sepa- 
ration might cause an azimuthal asymmetry about the shower axis labeled 
by a dipole p, then F(q) depends on F(q ■ p). The orientation of the dipole 
relative to the observation point thus has a strong effect on the emission. 
(Other numerically large effects will also be important: for example Allan|l]] 
incorrectly assumes the fields go like 1/R from his use of Feynman's formula) . 

As a consequence of separating out the form factor, the integrations have 
become effectively one-dimensional. 

3 Numerical Work 

At this stage we have a formula for the vector potential which is a product of 
a form factor and an object I ff (t],6) containing the information about the 
shower history. We will denote I FF (r),6) the Fresnel-Fraunhofer integral be- 
cause it interpolates between these regimes. In the Fraunhofer approximation 
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it is easily shown that the factorization is an exact kinematic feature of trans- 
lational symmetry as exemplified in Eq. (Q). If one makes a one-dimensional 
approximation, the Fraunhofer integral then evaluates the Fourier transform 



of the current 25, 26 



The factorization in the Fresnel zone is more demanding, yet should be 
an excellent approximation. When calculating I FF (r], 8), we cannot (as men- 
tioned earlier) consistently expand in powers of z'/R because 77 ~ 1 will be 
needed. 

It makes sense at this point to make a numerical comparison with previous 
work. Summarizing the results of an extensive Monte Carlo calculation in 
the Fraunhofer approximation, ZHS gave a numerical fit to the electric field 



R\E? HS (6)\ 7 v E r ,0-0 c , air V . 

with vq = 500MHz. This is the result of a global fit to many angles, energies 
and frequencies v = uj/2it < u . In this convention uo is positive. The 
normalized form factor is Fzhs(v) = 1/[1 + 0.4(i//z/ ) 2 ], as discussed below. 
In making the calculation, results for the field were also rescaled due to 
computer limitations. As a result, the field reported is strictly linear in the 
primary energy Eq. We will comment on this shortly. 

We calculated our own result proportional to I FF (rj = 0, 6) over a range 
of many frequencies and angles (Figure 2). Before doing the integral we scale 
out the electromagnetic and dimensional factors which are obvious. In our 
convention — 00 < uo < 00. The results of our numerical integration are quite 
well fit by 



R\E r,= ®{9) I e 1 

^— — = —aV2rrn max ujsm6exp [--(ka) 2 (cos6 - cos6» c ) 2 ] (10) 

r yjJ) C Z 

Putting in numerical values, this gives: 

R\EPr z °{6)\ 7 a n max v . 1 ,cos# — cos0 c . ,,, V , 

" 2.09 x 10- 7 --^-— -exp[--( 



F(lu) ' m 1000 GHz 1 2 V A(cos0) MHz 

(11) 

where 

, , 2 m GHz 
Afcos 6) = 0.048^ . 



>e a v 



( We have indicated that a is in units of m. and v in GHz.) Note the linear 
dependence on au, argued earlier to come from dimensional analysis applied 



11 



to the limit 77 — > 0. A cursory inspection shows that this result and the 
Monte Carlo have the same general features. 

To continue the numerical comparison we need numbers for the longitu- 
dinal spread parameter a and the number of charges at shower maximum 
n-max- There are several ways to estimate this. Running the ZHS code many 
times and fitting the output of a ITeV shower with a cutoff of 611KeV gives 
a = 1.5m, n max = 345. Using these and allowing for the factor of two in 
conventions gives agreement to a few percent in normalization with ZHS. 
However, the other way to do the calculation is to evaluate the product an max 
many times. This method is preferred because fluctuations in a and n max 
are correlated. Doing this gives an max = (570 ± 50)m at ITeV, which would 
predict a normalization factor of (1.2 ± 0.1) x 10~ 7 in Eq. (|S|). This (plus 
the angular dependence studied below) indicates that the factorized result is 
quite consistent with the Monte Carlo. 

In Figure 2 we show numerically integrated values of I FF Eq. (§). These 
factors appear directly in the fit just cited, and serve to check the formulas. 
The form factor has been divided out. For the range of parameters relevant 
to the problem, agreement is very good, and relative error is much less than 
1%. 

For experimental purposes one would like independent confirmation of 
the parameters from another source. Net particle evolution is well described 
by Greissen's classic solution, which was simplified further by Rossi to a 
Gaussian, n(z) = n max exp [— z 2 /(2a 2 )}/ \/2na. While Greissen's a refers to 
the whole shower, it should also be a reasonable description for the longi- 
tudinal spread of the charge imbalance, which tends to be a fixed fraction 
of the total number of particles after a few radiation lengths. (There is 
one caution that the Greissen formula does not explicitly include low-energy 
physics important for the charge imbalance.) The particular Greissen for- 
mula we consulted p7| for the longitudinal spread in radiation lengths X 
gives cl/Xq ~ ^j2/3log(E / E) for particles in the shower with energy greater 
than E and a primary with energy E Q . In that case one estimates a = 1.8m 
at Eq = ITeV with E = QllMeV, X = 0.39m in ice, which is acceptably 
close to the previous estimates. 

At higher energies there is every reason to believe that Greissen's 
stretched ^f\ogE energy dependence will apply. In that case a = 2.1m, 2.3m 
for lOOTeV, IPeV showers, respectively. Note that the product n max a is 
relevant for the field normalization. In this case we also need n max ~ 
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1/y log (Eq/E) — 0.33 from the same Greissen approximation. Rather amaz- 

ingly, the product n max a ~ (E /E) ^2/3 log {E /E) /^/\og {E /E) - 0.33 
~ Eq/E at high energies. This confirms the phenomenon observed by ZHS: 
the normalization of the electric field (Fraunhofer approximation, 6* = 9 C ) 
scales precisely linearly in the primary energy. It is rather pleasing that the 
result can be understood from first principles |21| . Later we will see that 
the parameter a enters in a much more complicated way in the Fresnel zone, 
creating an extra, weak energy dependence. 

Regarding the angular dependence, our work (Figures 3) indicates a gen- 
eral dependence on cosO — cos9 c rather than 9 — 9 C . When fitting numer- 
ical output the two functional forms are rather different, unless one has a 
very narrow distribution. Linearizing for small cos9 — cos9 c with a = 1.5m 
for the comparison, we would predict the scale in the angular dependence 
A9 = 2.1°{pq/p) while ZHS have the same expression with 2.4°. We find 
that A9 is proportional to 1/a. If a grows slowly with energy, as Greissen's 
formula indicates, then the angular width decreases, which is not seen in 
ZHS. Another possible explanation for the small discrepancy is the im- 
proper radiation from tracks terminating abruptly at the ends used in the 
Monte Carlo. We have identified these effects as responsible for the small 
oscillations seen in the Monte Carlo output, an effect apparently too small 
to measure. 

When numerical output to the frequency dependence is fit, there is a 
slight coupling between the model for the form factor and the parameter a 
one will extract for the longitudinal spread. We made our own fit to the 
ZHS code's frequency dependence including the region up to 1GHz, using 
a Gaussian form factor because of its better analytic properties. (The ZHS 
fit, which contains poles in both the upper and lower half-plane, violates 
causality.) Specifically, we find 

F[y) = exp [-v 2 /(2vl)]i = 0.93GHz. 

Using the corresponding a value we would predict the Fraunhofer Monte 
Carlo A9 ~ 2.3°(500MHz/v), quite close to ZHS. 

In real life, shower to shower fluctuations are highly important. We stud- 
ied the statistical features^ of the parameters n max , a by fitting individual 
showers many times and looking at the average and rms fit values. The re- 
sults at E — ITeV were a = 1.5 ± 0.2m, n max = 345 ± 60. Multiplying these 

2 We thank Soeb Razzaque for help with this. 
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and adding fluctuations in quadratures gives an max = 520(1 ± 0.22)m. The 
combination an max *, which is the primary variable in determining the nor- 
malization of the electric field, was found to be (570±50)m. The fluctuation 
of an max is less than half the value that the uncorrelated fluctuations of the 
separate terms would give. The relative fluctuation is said to decrease with 



increasing energy f28|j, but there are uncertainties. For example, threshold 
rescaling is used in Monte Carlos, leading to loss of information about the 
true fluctuations. Very preliminary results of running the standard Monte 
Carlo GEANT show variations in average shower parameters such as an max 
at the 30% level compared to the average of ZHS |JT]. These comparisons in- 



dicate that the electrodynamics is probably determined better than the rest 
of the problem. Indeed, the deviations from Gaussian behavior in showers 
is an effect contributing to the fields at the few percent level. At the level 
of 10 — 15%, many other small effects contribute. Unless one uses details 
about the uncertainties and errors in fits, and especially about the shower- 
to shower fluctuations in all relevant quantities, it is pointless to fine-tune 
the comparison further. We conclude from the numerical work that the fac- 
torized expression is at least as reliable as the Monte Carlos, and has the 
attractive feature that the parameters can be adjusted directly. 



4 The Saddle Point Approximation 

While the Rossi-Greissen Gaussian approximation to the shower is common, 
there are additional features which favor such an approach to the emitted 
radiation. The coherence is dominated by regions where the phases add 
constructively, greatly enhancing the peak region. In such circumstances 
analysis is helpful, especially when the largest contribution to I FF is dom- 
inated by saddle points. These are points where the phase is stationary, 
d(j)(z*)/dz* = 0. 

Here we describe the saddle-point method to evaluate I FF (r),9). This is 
a classic, controlled approximation when the charge distribution has a single 
maximum and kR ^> 1. The method turns out to give the exact result in 
the limit of flat charge evolution, that is, the Frank- Tamm formula, where 
numerical evaluation is highly unstable. With the saddle-point approxima- 
tion, we can extract analytic formulas which are as good as the numerical 
integration. 

We now describe the saddle-point features. By translational symmetry 
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the shower maximum can be located at z' = 0. Referring to the formula of 
Eq. (U) cited earlier, the saddle point of the phase is given by solving 

cos^c — (z — z*)/R(z*) — (i/k)(d/dz*) lognfz*) = 

for the point z* dominating the integral. The maximum electric field is 
already known to occur at points (z, R) near the Cherenkov cone. For such 
observation points z = Rcos9 c , and the saddle point equation has an easy 
solution at z* = 0. Thus the dominant integration region is near the shower 
maximum, as physically expected. 

As the point of observation moves off the Cherenkov cone, the saddle point 
moves away into the complex z' plane. To find the complex saddle-point, we 
approximate logn(z') ~ — z' 2 /(2a 2 ) in the vicinity of the shower maximum; 
that is, we fit the top of the shower locally with a Gaussian. To re-iterate: the 
saddle-point approximation does not need to replace the entire shower by a 
Gaussian, but replaces the vicinity of the region where phases are contribut- 
ing coherently by a Gaussian. The saddle-point condition gives a quartic 
equation which can be solved. Unfortunately the solution is impossibly com- 
plicated, thwarting a direct approach. We circumvented this by studying 
the saddle-point location numerically. We found the quartic solution is accu- 
rately linearized in a special variable: expand about cos# close to cos# c . We 
then find z„ « R(0) sin 2 fl(cos 9 - cosfl c )[l + iR(0)/ka 2 sin 2 9}-\ With this 
formula the reader inclined can repeat all the calculations. We also show the 
saddle-point to highlight the appearance of the ratio R(0)/(ka 2 sin 2 9) indi- 
cated by the qualitative "acceleration" argument. Finally, replace R(0) = 
R = Vz 2 + p 2 . 

The rest of the calculation is standard mathematical physics [p9| |, so we 
just quote the results. Calculating fields from the potential and keeping only 
the leading terms in l/(kR) C 1 we find: 

K = ^F(q)I FF ( V ,9)[(cos9 -cos9 c )e R 

. cos 9 C cos 9 - cos 9 C . 
-t 1 -' 17 !—^ 1 : )sm0ee], (12) 

5 iuJ vf^rFF, fl vi i • COS ^ COS#-COS# c . 

B " = 5 T F ^ 1 (V,0)(1 + IV— — a ; : )smfle , (13) 

vcR cos 9 C stn z 9 1 — trj 

where 

tPI7 , „. ;i . R / — . . n . . cos 9 cos 9 — cos 9 C 1/9 

I FF ( V , 9) = e lkR aV2^[l - i V (8)(l - 3t V -^- : c - )]- 1 ' 2 

sin 9 1 — ir] 
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. 1,, x9 (cos 9 — cos6U 2 -, . . 

exp[--(fai) 2 * — ], (14) 

and where 

/ca .2/1 
77 = — — sm 

Inspection reveals that these formulas have the dependence on u, R, 
and a distance scale a quoted earlier on physical grounds. The formulas 
cited earlier as summarizing the numerical work are, in fact, the saddle- 
point approximations evaluated at 77 = 0, which fit closer than any empirical 
formula. 

4.1 Remarks 

The dependence of the fields on symbol 77 summarizes a good deal of com- 
plexity. For example: 

• The limit rj — > yields the Fraunhofer limit, with spherical wave fronts 
and Eu ~ aui / R. 

• The limit 77 — > 00 gives the cylindrically symmetric E u ~ J~lo / R fields. 
This field can be substantially different from the Fraunhofer approxima- 
tion. In fact one must take this limit to get the Frank- Tamm formula. 

A notable application is emission from ultra-high energy air showers. 
In such showers the LPM effect plays a definite role in suppressing 
the soft emissions from the hardest charges. There is no corresponding 
suppression of the evolution of the low-energy regions of showers where 
most particles exist, however JjlJ, [18| . The major effect that we find is 
that the showers become long kinematically: that is, the parameter a 
gets big if one is working with, say, a primary of 10 20 eV. The emitted 
fields approach those of the Frank- Tamm formula as a — > 00, via the 
formula cited. The effect is important numerically: a 10 20 eV air shower 
does not approach the Fraunhofer limit nearer than 300km. The con- 
ditions of RICE are more amenable to the limit, and at R ~ lkm, 
lo ~ 1GHz the Fraunhofer-based estimates near the Cherenkov cone in 
the literature are good to 20%. 

• The frequency dependence in the Fraunhofer approximation is strongly 
affected by "diffraction". Independent of the form factor effect, the 
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Fraunhofer approximation imposes an upper limit to the frequency of 
order uj < (1/ a) (cos 6 — cos# c ). The true behavior is substantial 
different: from Eq. (|I~4| ) we see that the field exists in a region 

uj < (l/a)[(cos0 - cos9 c ) 2 - (a 2 /R 2 ) sin 4 (#)]- 1/2 . 

For large u and in the angular region where the signal exists, the be- 
haviour is much flatter. This is illustrated in Figures 4 and 5. This 
effect is invisible at the exact value 9 = 9 C , where the fake Fraunhofer 
frequency cutoff and most of the true functional dependence in the ex- 
ponent of I FF both drop out. As a result of the difference in frequency 
dependence, the time-structure of the electric field may be substan- 
tially different from the Fraunhofer approximation. We will return to 
this point in a Section Causal Features . 

Fields in the forward and backward directions, sin 2 # — » 0, are the fields 
of f] — > 0, the Fraunhofer approximation, regardless of the physical val- 
ues of k, a, R. We note that the experiments of Takahashi [[13] observe 
an extremely limited region of sin 2 # — > 0. Perhaps this contributes to 
the observed agreement with Tamm's formula in a regime where a 2 k/R 
is not close to zero. 

The polarization varies considerably. From symmetry the polarization 
is in the plane of the charge and the observation point. Moreover, 
for 9 = 9 C the electric field is transverse to R for any rj. Yet naive 
transversality is not true in general at any finite rj. 

Between the various limits the dependence on every scale in the prob- 
lem, namely the frequency u, the distance R, the length scale a, and the 
angle 9, is neither that of the Fraunhofer limit nor that of the infinite 
track Frank- Tamm limit, but instead a smooth interpolation between 
the cylindrical and spherical wave regimes. 

In the finite rj limit, one may also include a further effect, namely 
that as uj — > one has a 'near-zone' Coulomb-like response at small 
R. (Indeed, the uj — > limit measures the net charge.) This effect, 
important below about 10MHz, also has a slight effect on the time 
structure of pulses. 
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We pause to comment on the generality of the result. What if we had 
not made the physical, but specific ansatz (j3|)? The entire analysis can be 
repeated for an arbitrary charge distribution j(t', z', p'). The Fraunhofer ex- 
pansion of the transverse variable, and the Fourier integral of the t' variable 
are general. Provided the p' extent is finite, and there exists a dominant z' 
region, then the integrals always factor into a product of a form factor and 
a one-dimensional integral for A^ F (rj). In fact nothing changes (the reader 
can repeat the calculation) except that when an arbitrary current is set up, 
the existence of a single saddle point cannot be assumed. Corrections to the 
local Gaussian approximation are straightforward. The slight skewness of 
real showers (or other arbitrary charge distributions) can also be developed 
as a saddle-point power series. Again: there are elements of bremmstrahlung 
in real showers, having a stochastic nature, which the current model has not 
attempted to reproduce. Detailed Monte Carlo simulations |JT| of our group 



have included the bend-by-bend amplitudes of tracks undergoing collisions. 
This goes well beyond the approximation of a single straight line track, sud- 
denly beginning and ending, of the previous literature [jT^, jlBfl . The effect 
of all the small kinks is negligible except in the very high frequency region 
u) ^> lOOG-f/z, while the endpoint accelerations give oscillations in the an- 
gular dependence down by orders of magnitude. As a final side remark: 
we explicitly studied contributions of finite tracks, just to see what would 
happen, in the development towards the conditions of the Tamm formula. 
It is straightforward to develop these pieces if one needs them for, say, the 
Takahashi-type experiments El, in the Fresnel zone. 



4.2 The General Case 

We now turn to fields valid for any r/. The form factor, which was extracted 
from the Fraunhofer calculations, is universal and need not be changed. The 
validity of the saddle-point approximation does not depend explicitly on the 
value of r\. The procedure of linearization to locate the saddle-point happens 
to be good to cosO — cos9 c ~ 1, so that the approximation is rather good in 
the entire region i?/a> 1, kR 3> 1. 

For practical applications it is useful to have a formula for the fields 
with quantities measured in physically motivated units. For this purpose we 
rewrite flT2|) as follows 

* 4K2 . 52 .nn^^( 9 -)^[^]. d5) 
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Here n max is the excess of electrons over positrons at the shower maximum, 
R is measured in meters, v is measured in GHz. The rescaled field is 

rf cos# — cos# c ^ , . cos 9 r cos 9 — cos 9 r N 
£ = [ — n £r + (1 - ir]- 



sin 9 " sin 2 9 1 — %r\ 

cos 9 cos 9 — cos 9 C . . , /9 

1 _ ^ 1 - 3^-^- ; ^ - 7 

sur9 I -it] 

r !/, x 2 ( cos ^ - cos^ c ) 2 
exp[--(fca) ^— — ]. (16) 

We have defined a kinematic factor ^ = —iexp(ikR) sin 9 in such a way that 
the rescaled field 8 is normalized at 9 — 9 C , 

8(9 = 9 c ) = (l-i7])- 1 / 2 eg. (17) 

It is convenient to plot the magnitude of the rescaled field, Eq. flTBT). 
Figure 3 shows the magnitude as a function of the angle difference 9 — 9 C in 
various limits. The Fraunhofer approximation is shown by a dashed curve, 
and our result by the solid curve. One observes that the Fraunhofer limit is 
approached from below. This is physically clear: The Fresnel zone fields have 
a wider angular spread, and conservation of energy forces them to be smaller 
in magnitude compared to the sharper, diffraction-limited Fraunhofer fields. 
As the fields evolve to infinity, they coalese into narrower and taller beams. 

The frequency dependence of R\E W \/ F{ui) is shown on Figures 4 and 
5. Exactly at the Cherenkov angle the difference between the Fraunhofer 
approximation and our results are minor for the typical parameters of RICE. 
However, away from the Cherenkov angle there is a substantial difference 
between the two, throughout the region where the magnitude of the field 
is large. This effect can be masked by the form factor, so we have plotted 
RE^j ' F{lS) to show it. This effect may have important repurcussions for the 
time-structure of pulses, which are also discussed in the last Section. 

Figures 6-10 are contour plots of the electric field. We did not bother 
to remove the small region a/R ~ 1, where our result does not apply. The 
Fraunhofer approximation has trivial 1/R dependence on the distance to the 
observation point (Figure 10). The exact result is certainly different, with 
Figures 6, 8, and 9, in particular, illlustrating the effects of constructive 
interference in the region of cylindrical symmetry. A complementary view 
examines contour plots of constant phase. This is shown in Figures 8 and 
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9. In making these figures we decreased the kinematic phase ikR to val- 
ues showing several oscillations (as opposed to hundreds) across the range 
of the plots. The lines of constant phase can be used to illustrate the time- 
evolution of waves of a given frequency: that is, the Fourier transform of 
5{uj — uj*)Euexp{—iujt) has wave fronts at each moment in time given by the 
lines of constant phase. The constant phase lines of the Fraunhofer approxi- 
mation are, of course, spherical (Figure 10). The constant phase lines of the 
true behaviour interpolate between cylindrical and spherical symmetry. As 
a consequence of the Fresnel-zone behaviour, the eikonals of the expanding 
radiation field do not emerge radially, but actually curve due to interference 
effects. This is a sobering impact of very basic physics, which has a measur- 
able effect in the signal propagation speed discussed later under the topic of 
causal feature. 

5 Causal Features 

With our convention that the electric field E(t, x) = J^dto exp(-iujt) 

E u (x), causality requires E u to be analytic in the upper half-plane of u. 
Singularities in the lower half-plane determine the details of E(t,x) and the 
precise causal structure. 

To discuss this we consider detection of signals via an antenna-system 
response function A w . By standard arguments the detected voltage is a 
convolution in time, and therefore a product in uj space, of the antenna 
function and the perturbing electric field. The antenna function has the 
same causal analytic properties as the electric field. A proto-typical antenna 
or circuit function for a driven LRC circuit is 

Au = Z/{-u 2 + uj 2 a - iTu), 

where Z depends on where the amplifier is connected in the circuit and can be 
treated as a constant. Note that A u = —Zj (cu — ou + )(cu — uJ) where uj± are in 
the lower half-plane. The dielectric function e(u) can then be taken as slowly 
varying in the region where the antenna and form factor allow a response, 
and also has its analytic structure in the lower half-plane if this detail needs 
to be included. We will also ignore the form factor for this discussion, which 
earlier was cited as a formula analytic in the complex plane. While nothing 
in our analysis depends on these idealizations, this approach to the analytic 
structure serves to make our point. 
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As a first illustration, consider the electric field fit given by ZHS, pro- 
portional to — iuJo){uj + iu Q ) with u /(2it) = 500MHz. This field has 
poles in both the upper and lower half-plane, violating causality. One may 
argue that the literal analytic behavior in the complex plane goes beyond the 
ambitions of the original semi-empirical fit. Nevertheless, Cauchy's theorem 
applies to the subsequent numerical integrations that have been made ||, 
giving non-causal branches to numerically evaluated Fourier transforms, as 
well as unphysical short-time structure. 

Let us compare the analytic structure of the electric field in the saddle- 
point approximation. This approximation does not attempt to describe the 
region u — > 0, which requires treatment of the near zone. However for 
causality we do not need E u near the origin but at large |a>|. The saddle 
point approximation is good here so the results should be reliable. 

Let us investigate this in more detail. In the exponent in the expression 
for the electric field we have — (1 / 2) (ka) 2 (cos 9— cos 9 C ) 2 (1+ir/) / (1+r/ 2 ) . Since 
rj goes like uj there is a phase linear in uj at large |a>|. There is also a branch- 
cut and pole from the prefactor which occurs at rj — — i. All singularities are 
in the lower half-plane and consistent with causality. The causal structure 
of E(t,x) then hinges on closing the contour at infinity. For this the details 
of the antenna function, which generally has isolated singularities, as well as 
singular behavior of E u near the origin do not matter. 

We close the contour at infinity avoiding the branch cuts, which may be 
oriented along the negative imaginary axis or as consistent. Convergence 
then requires 



This result has a natural interpretation. While the radiation from the 
shower appears to come primarily from the geometric location of the maxi- 
mum, 9 ~ 9 C , the shower actually develops and radiates earlier. Consequently 
the strict causal limit must correspond to an apparent propagation speed 
slightly faster than the naive speed of c/ yfe deduced from the location of the 
observation point at R. The earliest signal actually arrives at an apparent 



lim Re[— icut + ituy/eR o> 2 ea 2 (cos# — cos# c ) 2 -] < 0. 

*oo 2 1 + r/ z 



Using the definition of 77, causality implies 
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speed v app of 

c (cos9 — cos9 c ) 2 ! 
y/e 2 sin 6* 

This formula is entirely geometrical, consistent with the simple picture that 
distance differences in the problem are causing the effect, but it also incorpo- 
rates subtle features of coherence. For example, the distance scale a cancels 
out. Yet a determines the angular spread cos 9 — cos 6^ over which most of 
the power in the wave is contained, and enters in this fashion. 

The Fraunhofer limit, in which all signals originate at a single point of 
origin, is incapable of capturing such an effect. It is interesting to trace the 
origin of the discrepancy. The singularities of interest are located by the 
zeroes of 1 — ir\ = 1 — iu \fe sin 2 9 / R. When the limit R — > oo is taken in 
the first step of the Fraunhofer approximation, all the non-trivial analytic 
structure moves away to uo — > — too and is lost. This procedure does not 
commute with closing the contour at — > oo. The correct procedure, of 
course, is to first close the contour, and then take the limit of large R. 

In practice, of course, not all of the signal arrives at the earliest possi- 
ble moment. The time scale over which the signal is detected depends on 
competition between dispersion, the antenna and form factor details, and 
something like twice the "advanced" time interval. This time interval is 
At caus = i?y / e/(2c)(cos6 l — cos9 c ) 2 / sin 2 9. Since the At caus effect scales pro- 
portional to the distance R, it does not become negligible in any limit, ex- 
hibiting another subtle facet of breakdown of the Fraunhofer approximation. 
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Figure 1: Cartoon of electromagnetic coherence zone associated with an 
evolving shower. The shower develops over a longitudinal spread dimension 
a, as represented by the smooth curve. During most of this development the 
pancake of charge has a constant size d, but an evolving charge normaliza- 
tion, represented by the gray scale. The coherence zone Az co h depends on 
the angle, frequency, and location of the observation point. The situation 
illustrated has Az co h > a, a case in which full coherence of the charge (up 
to a frequency cutoff determined by the pancake size) is obtained. In the 
reverse limit Az co h < a, only the portion of the charge inside the region of 
Az co h contributes constructively. 
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Figure 2: Comparison of numerically integrated (points) and analytic fit 
(solid curve) to the Fresnel-Fraunhofer integral I FF (r],9), Eq. @. Values 
of parameters are: a = 1.5m, R = 1000m; frequency v = 1,2,5,10GHz 
(from top to bottom). Both Fresnel and Fraunhofer regimes are successfully 
reproduced. Angle is measured in degrees, I FF in meters. 
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Figure 3: Magnitude of the rescaled electric field E = \£\, Eq. Our 
result (solid curve) is compared to the Fraunhofer approximation (dashed 
curve), a = 1.5m, R = 50m, v = 1GHz; angle is measured in degrees. 
The Fraunhofer approximation is narrower in angular extent and larger in 
magnitude. Due to conservation of energy the field sharpens its angular 
distribution and approaches the asymptotically far field from below as R — > 
oo. Up to scaling normalization factors, the same plot would apply to any 
with the same m for example a = 3m, R = 200m. 
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Figure 4: Magnitude of the re-scaled electric field W = 
(1000/n max )R\E UJ \/F( y uj) as a function of frequency u, evaluated at 
6 = 9 C + 1°. W is measured in 10~ 7 V/ MHz. The dashed curve is the 
Fraunhofer approximation. Note the high frequency cutoff, which is not a 
form factor effect, but an artificial result of the limit R — > oo imposed by 
the approximation. Solid curve is our result. The parameters a = 1.5m, 
R = 100m. 
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Figure 5: Magnitude of the re-scaled electric field W = 
(1000/n max )R\E UJ \/F(uj) as a function of frequency is, evaluated at 
6 = 9 C + 0.3°. W is measured in 10~ 7 V/MHz. The dashed curve is the 
Fraunhofer approximation; solid curve our result. As R — > oo the Fraunhofer 
approximation begins to apply. The parameters a = 1.5m, R = 1000m. 
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Figure 6: Contour plot for the magnitude of the electric field, a = 5m, 
v = 100MHz, distances in m. Note the evolution of the field from cylindrical 
to spherical behaviour as the distance from the origin increases. 
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Figure 7: Contour plot for the magnitude of the electric field in the Fraun- 
hofer approximation. The parameters a = 5m, v = 100MHz; distances in 
m. The magnitude of the field lacks the richness of structure of Figure 6. 
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Figure 8: Contour plot for the phase of the electric field. The parameters 
a = 1.5m, v = 100MHz; distances in m. Each contour represents the time 
evolution of an expanding wave front for the frequency used. The phases 
illustrate non-trivial propagation with destructive interference, as well as 
evolution from cylindrical to spherical symmetry. 
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Figure 9: Contour plot for the phase of the electric field. The parameters 
a = 1.5m, v = 100MHz; distances in m. The phases illustrate non-trivial 
propagation with destructive interference, here shown on a larger scale than 
Figure 8. 
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Figure 10: Contour plot for the phase of the electric field in the Fraunhofer 
approximation with a = 1.5m, v = 100MHz; distances in m. The contours 
illustrate kinematic spherical symmetry imposed by the approximation. 
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